Short range Coulomb correlations render massive Dirac fermions massless 
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Tight binding electrons on a honeycomb lattice are described by an effective Dirac theory at low energies. 
Lowering symmetry by an alternate ionic potential (A) generates a single-particle gap in the spectrum. We 
employ the dynamical mean field theory (DMFT) technique, to study the effect of on-site electron correlation 
(U) on massive Dirac fermions. For a fixed mass parameter A, we find that beyond a critical value Uci{A) 
massive Dirac fermions become massless. Further increasing U beyond (7c2(A), there will be another phase 
transition to the Mott insulating state. Therefore the competition between the single-particle gap parameter. A, 
and the Hubbard U restores the semi-metallic nature of the parent Hamiltonian. The width of the intermediate 
semi-metallic regime shrinks by increasing the ionic potential. However, at small values of A, there is a wide 
interval of U values for which the system remains semi-metal. 

PACS numbers: 73.22.Pr, 71.30.+h, 
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INTRODUCTION 

The single-particle spectrum of excitations in graphene 
can be very accurately described by a tight-binding model in- 
volving hopping between the localized orbitals of neigh- 
boring carbon atoms The low-energy sector of such 
Hamiltonian will be described by a 2h-1 dimensional Dirac 
theory 130. Starting from this point, many perturbations can 
be imagined and/or fabricated to modify the spectrum of Dirac 
electrons From applications point of view, it is important 
to open up a gap in the spectrum by lowering the symmetry 
of the nearest neighbor tight-binding Hamiltonian. For ex- 
ample the substrate can induce a sub-lattice symmetry break- 
ing, e.g. by adding an alternate ionic potential of strength A, 
which immediately leads to a charge gap of magnitude 2A 
in the spectrum of single-particle excitations and renders the 
Dirac electrons massive |5, 62 the ground state of which will 
be a trivial band insulator (BI). Recent ab-initio estimates of 
the strength of the Hubbard U in graphene suggests that the 
on-site Coulomb repulsion is quite remarkable, ^ 10 eV 101. 
Since this value is expected to be a local property of pz or- 
bitals in carbon atoms, one do not expect the above value of 
U to be much different when the Dirac fermions of the un- 
derlying honeycomb lattice acquire a mass due to extrinsic 
effects, such as substrate, binding with ad-atoms such as hy- 
drogen ISO, etc. Therefore it is important to consider the com- 
petition between the single particle gap parameter A and the 
Hubbard parameter U on top of a semi-metallic state at half- 
filling. On the strongly correlated limit, when the Hubbard 
energy scale dominates over the ionic potential, i.e. U ^ A, 
the system will be a Mott insulator (MI) where the charge fluc- 
tuations become costly because of no-double-occupancy con- 
straint imposed by large U, while in the opposite limit where 
A dominates over U, the system can be described in terms of 
an effective massive Dirac theory. The purpose of this paper 
is to show show that short range many-body Coulomb inter- 
action U can transform massive Dirac fermions into massless 



ones. 

The threefold coordination of carbon atoms on the hon- 
eycomb lattice of graphene is the basic mathematical reason 
for emergence of cone-like dispersion in graphene which is 
responsible for semi-metallic properties. If the dimension 
and/or coordination number were different, the parent one- 
band tight-binding Hamiltonian would at half-filling describe 
a metal. Previous DMFT study by Garg and coworkers pre- 
dicts that for metallic parent Hamiltonians, the intermediate 
phase will be a metal [9J, i.e. the competition between Hub- 
bard and ionic terms in a metal, restores the metallic state. 
One may ask the same question for semi-metallic (Dirac) sys- 
tems; What would be the result of competition between the 
Hubbard and ionic terms in a semi-metal? Would the inter- 
mediate phase be a metal? Or the competition between the 
Hubbard and ionic terms -both of which would individually 
drive the system towards insulating behavior - will give way 
for the semi-metallic behavior? 

To tackle this question, the dynamical mean field theory 
(DMFT) is a suitable and powerful method fl^. which is 
capable of handling Hubbard interaction in a proper way. 
This approximate technique becomes exact in the limit of 
infinite coordination numbers ifiol, [Till which can be for- 
mally incorporated by assumption of a higher dimensional 
generalization of various lattices [12I1 . For lower coor- 
dination number, the local (k-independent) self-energy em- 
ployed in the DMFT becomes only an approximate descrip- 
tion. Hence, the most significant drawback of the method 
is expected to be the underestimation of the spatial quantum 
fluctuations Hs- 151. Therefore the values of the critical pa- 
rameters obtained from a simple DMFT on the honeycomb 
geometry maybe overestimated US\ . However, the overall 
picture emerging from this numerically powerful method is 
expected to hold even when more sophisticated methods, such 
as cluster DMFT 1I17I1 . Gutzwiller projection llSll alterna- 
tive techniques such as quantum Monte Carlo lll9ll or slave- 
particle |2QI1 schemes are employed. Since in this work we are 
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not interested in fine spectral details, the lack of k-resolution 
we believe will not affect the central conclusion of the paper, 
namely that the semi-metallic behavior is restored when Hub- 
bard and ionic potential compete. 

The ionic-Hubbard model has been studied extensively in 
one dimension L21i1 . and two dimensions, and the nature of 



matrix form as. 



interim phase has been debated Il22ll . In addition to various 
low-dimensional techniques employed to study such model, 
the DMFT technique has also been employed to study the na- 
ture of intermediate phase of this model for metallic parent 
systems flEl. In Ref. " 



19(1 the authors implemented DMFT 
for the Bethe lattice and found that by increasing U for fixed 
A, the system undergoes two phase transition: First from band 
insulator into metallic phase, and second from metal into Mott 
insulator phase |.'>|. However in Ref. ll23ll they used almost the 
same method and find coexistence insulator phase region for 
some A,U 



MODEL AND METHOD 

The ionic-Hubbard Hamiltonian on the honeycomb lattice 
with two atoms in same unit cell is given by. 



H=-t [cl^cj,^ + h.c] + 

i£A,j£B 



(1) 



where t is the nearest neighbor hopping, A denotes the ionic 
potential which alternates sign between site in sub-lattice A or 
B. The Hubbard U stands for the repulsion potential between 
two electrons with opposite spins in ith site of the lattice. The 
chemical potential is /i = C//2 at half filling and so the aver- 
age filling factor is — \ xhe model for t > Q and 
noninteracting limit [/ = 0, represents a band insulator, with 
energy gap Egap = 2A. The average filling factors at half- 
filling will become (ua) = 0, {ns) = 2. For U ^ A, the 
system is in the Mott insulator phase with ua = = 1 124 1. 



In the intermediate region, a flow-equation analysis shows that 
at some energy scale, the ionic and Hubbard terms are ex- 
pected to cancel each other's effect 12411 . From this angle we 
expect the character of the parent Hamiltonian (i.e. the hop- 
ping term only) to dominate the nature of the ground state. 
Therefore for metallic parent Hamiltonians we expect the in- 
termediate phase to be a metal ||9, 24], while for semi-metallic 
parent Hamiltonian as will be shown in the following, the 
semi-metallic character will be restored and the massive Dirac 
fermions will become massless. as a result of increasing the 
Hubbard U. 

To formulate the DMFT machinery for semi-metallic par- 
ent system, consider a parent tight-binding Hamiltonian on 
the honeycomb lattice to be in the paramagnetic phase. The 
interaction Green's function in the bipartite lattice acquires a 



G(fc,a;+) 
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where k is the momentum vector in first Brillouin zone, e(fc) 
is the energy dispersion for the honeycomb lattice 13], and 
Ca(B) = T A + M ^ ^A{B){^^) with CJ+ = w + zO+. In 
DMFT approximation the self-energy, T.a{uj+) is local ifioll . 
i.e. the self-energy matrix is diagonal and independent of k. 
Therefore the off-diagonal elements vanish. Hence the local 
Green's function corresponding to sub-lattice a = A,B can 



be written as. 
tola, 
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J2k^aa{k,i-o'^) which simplifies 



de 



Po(e) 



CA{uJ+)CBiuJ+) - 



(3) 



where a = A{B), a = B{A) and pQ{e) is the bare DOS of 
the honeycomb lattice Q]- 

We start with an initial guess for the filling factor and self- 
energy Then we determine host Green's function from 
the Dyson's equation 5t7j(a;+) = G~^(a;+) + Y,a{uj+). Af- 
terwards, we solve impurity problem and find S]Q(a;+) — 
^aiGoaii^^)]- In this step we use iterated perturbation the- 
ory (IPT) as impurity solver ifioll . The iteration of these 
steps continues until convergence is reached. After conver- 
gence, we calculate the density of states given by Pq(w) = 
— Im Tr[GQ,(fc, w+)]/7r. The particle-hole symmetry at 
half-filling leads to pa{^) ~ Pb{~^), for the DOS of the 
two sub-lattices. The total DOS for the honeycomb lattice is 
eventually obtained via p{uj) = pa{'-^) + Pb{^)- 



RESULTS AND DISCUSSION 
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FIG. 1 : (Color online) The density of state for honeycomb lattice in 
ionic-Hubbard model is plotted for U = O.lt, It, lit and A = Q.lt 
from down to up. Between band and Mott insulator phase the semi- 
metal phase exist as the middle phase. 

In Fig. [T] we have shown the converged results for three 
typical values ofU at a constant A — Q.lt. As can be seen in 
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FIG. 2: (Color online) The zoom in the low energy part of the DOS 
for U = O.lt, It, lit at fixed A = O.lt. By increasing U , the single 
particle gap closes to restore the semi-metallic phase, which is then 
followed by a SMIT. 



this overall picture which covers the whole range of energies 
in the bandwidth, for small values of U, we have simple band 
gap in the spectrum which has single-particle character As 
U increases, spectral weight is transferred to higher energies 
to form upper and lower Hubbard bands. The striking feature 
in Fig. [T] is that, despite the formation of Hubbard bands at 
larger values of U, the overall shape of the low-energy spectral 
features resemble that of the parent Hamiltonian. The DOS 
around the Fermi level is V shaped at energy scales above 
^ A, and the van-Hov singularity arising due to the saddle 
point in the band structure moves towards lower energies. 

In Fig. 12] we have zoomed in the lowest energy scales 
around the Fermi level. As can be clearly seen in this fig- 
ure, increasing U eventually closes the single-particle gap, 
and restores the linear DOS at all low-energy scales, including 
those below ^ A, and renders the system semi-metallic. To 
our knowledge, this is a rare instant where many-body inter- 
actions of relativistic fermions ejfectively removes their mass 
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FIG. 3: (Color online) By increasing U the slope of DOS in the 
intermediate SM region changes. This in turn means an increase in 
the Fermi velocity of electrons. This figure corresponds to A = G.lt. 



and makes them massless. By further increasing U, up to lit, 
the spectral transfer to upper and lower Hubbard bands be- 
comes more and more (Fig.[Tli, while at the same time a Mott- 
Hubbard gap appears around the Fermi surface as in Fig. |2] 
Therefore at a fixed value of A, as one increases U, there will 
be a phase transition at Uci to the parent semi-metallic phase, 
followed by another phase transition at larger Uc2 to the Mott 
insulating phase. The intermediate phase is characterized by 
a strictly linear DOS in the low-energy scales, which quali- 
fies the intermediate phase for an effective Dirac theory. Note 
that the Fermi velocity in the intermediate effective 2h- 1 Dirac 
theory is renormalized with respect to the effective Dirac the- 
ory of the parent Hamiltonian. In Fig. [3] we show converged 
results for three typical U values in the intermediate regime 
which indicate the decrease in the Fermi velocity of the in- 
termediate Dirac liquid phase as a function of U which ulti- 
mately leads to a SMIT 111. 

Repeating such calculations for a range of ionic potential 
A, we map out the phase diagram for the ionic Hubbard model 
on the honeycomb lattice, the boundaries of which are repre- 
sented by blue and red lines in Fig.|4] First of all for A « 0, 
the width of the band insulating (BI) phase diminishes which 
is equivalent to Uci ~ 0, and there will be only one phase tran- 
sition separating SM and MI phases |[T3-15]. Note that the 
vertical axis denotes values of U /W, where W is the band- 
width of the parent Hamiltonian. As can be seen in Fig.|4] by 
increasing A, the SM region becomes narrower, and beyond 
A/i ^ 0.5 - the energy scale up to which the Dirac cone lin- 
earization holds - the red and blue boundaries merge. Beyond 
this point, the phase transition will be essentially between the 
BI and MI phases. 

For comparison in Fig.|4]we have also reproduced the data 
from Ref. 191] for a metallic parent Hamiltonian with Bethe lat- 
tice DOS. One interesting difference between the intermediate 
region on the honeycomb and Bethe lattice is that the in the 
former case, the ground state of the parent Hamiltonian (SM) 
occupies a much larger region in the two dimensional space 
of ( Y, ■^). This could be interpreted as the robustness of the 
parent phase of the Dirac electrons compared to Schroedinger 
electrons when an ionic and Hubbard energy scales compete. 
One more conclusion could be drawn from Fig.|4]is that if we 
imagine a horizontal line at a constant Ug/W ^ 0.5 — 0.7 
which maybe relevant to graphene |01, it can be seen that for 
a range of A/i ~ 0.04 — 0.05, the system still remains in 
the SM phase. For a reasonable estimate of t 2.7 — 3 eV, 
this implies that, graphene can remain SM despite a symmetry 
breaking ionic potential of strength A ~ 110 — 150 meV. Such 
a single-site DMFT may have overestimated the upper bound- 
ary C^c2(A) |jT3l. Improving on this calculations by e.g. clus- 
ter extension of DMFT is expected to push the upper boundary 
down, and hence our estimate of the tolerance A ^ 110 — 150 
meV is not expected to change much. 

Normally the gap magnitudes in gapped graphene samples 
are assumed to satisfy EqI^ = 2A. However this picture 
holds when one completely ignores the many-body effects. 
Hence we have used the superscript "0" to emphasize that it is 
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FIG. 4: (Color online) The DMFT phase diagram of the ionic- 
Hubbard model on the honeycomb lattice. Band insulating phase 
(BI) and Mott insulating phase (MI) are separated by a semi-metallic 
phase (SM). For comparison we have also given the similar phase 
diagram obtained for metallic parent Hamiltonian f^. The Hubbard 
energy scales in both plots are reported in units of bandwidth W. 



energy theory is 2+lD massive Dirac theory. For Uci{A) < 
U < C/c2(A) region the competition between U and A en- 
ergy scales restores the semi-metallic character of the parent 
Hamiltonian and makes them massless Dirac fermions. Ulti- 
mately, for U > [/c2(A) the strong correlations transform the 
intermediate SM to a Mott insulator lljl. Beyond A ~ 0.5t, 
the upper and lower critical values approach each other and 
the semi-metallic region will be so narrow that it will dis- 
appear This theoretical study has implications for gapped 
graphene samples. The observed values of spectral gap is al- 
ways below the "non-interacting" value of 2A. The strong re 
normalization of the gap magnitude is not the only manifesta- 
tion of the Hubbard correlations U. It also strongly influences 
other spectral features such as the life-time of quasi particles, 
e.g. in hydrogenated graphene OSll . 
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FIG. 5: (Color online) The Scap in units of t as a function of U/t 
for A = O.lt is plotted. Presence of non-zero U strongly modifies 
the actual gap, compared to the "non-interacting gap" 2A. 



"non-interacting" gap. In Fig.|5]we plot the actual spectral gap 
as a function of U for a fixed value of A = O.lt. As can be 
seen the actual (interacting) gap strongly depends on the inter- 
action parameter U in both band- and Mott-insulating phases. 
In particular the value of gap corresponding to U ^ At which 
is relevant to graphene, at A = O.lt ^ 280meV is expected 
to be less than half of the non-interacting estimate. Hence 
without proper account of the on-site correlations the es- 
timates of A based on the observation of -Ecap is expected to 
be about a factor of two underestimated. 



SUMMARY AND DISCUSSIONS 

We studied the ionic-Hubbard model on honeycomb lattice 
by DMFT method, and we found the system under goes two 
phase transition by increasing U for A > 0. For U < Uci{A) 
the system is in band insulator phase whose effective low- 



[1] A. K. Geim, K. S. Novoselov, Nature Materials, 6 (2007) 183. 
[2] A. H. Castro Neto, F. Guinea, N. M. R. Peres, K. S. Novoselvo, 

A. K. Geim, Rev. Mod. Phys. 81 (2009) 109. 
[3] G. W. Semenoff, Phys. Rev. Lett. 53 (1984) 2449. 
[4] M. Y. Han, B. Ozyilmaz, Y. Zhang, P Kim, Phys. Rev. Lett., 98 

(2007) 206805. 

[5] P K. Pyatkovskiy, J. Phys. Condens. Matter, 21 (2009) 025506. 
[6] G. W. Semenoff, V. Semenoff, Fei Zhou, Phys. Rev. Lett. 101 

(2008) 087204. 

[7] T. O. Wehling, E. Sasioglu, C. Friedrich, A. I. Lichtenstein, M. 

I. Katsnelson, S. Blugel, Phys. Rev. Lett., 106 (2011) 236805. 
[8] J. O. Sofo, A. S. Chaudhari, G. D. Barber, 75 (2007) 153401; 
[9] A. Garg, H. R. Krishnamurthy, M. Randeria, Phys. Rev. Lett. 

97, 046403 (2006). 
[10] A. Georges, et al. Rev. Mod. Phys. 68, 13 (1996). 
[11] M. Caffarel, W. Krauth, Phys. Rev. Lett. 72, 1545 (1994). 
[12] G. Santoro, M. Airoldi, S. Sorella, and E. Tosatti, Phys. Rev. B 

47(1993) 16216. 
[13] S. A . Jafari, Eur Phys. J. B. 68, 537 (2009). 
[14] M. Ebrahimkhas, Phys. Lett. A. 375, 3223 (2011). 
[15] M. T. Tran, K. Kuroki, Phys. Rev. B. 79, 125125 (2009). 
[16] Wei Wu, Yao-Hua Chen, Hong-Shuai Tao, Ning-Hua Tong, and 

Wu-Ming Liu, Phys. Rev. B 82 (2010) 245102. 
[17] A. Liebsch, Phys. Rev. B. 83 (2011) 035113. 
[18] L. M. Martelo, M. Dzierzawa, L. Siffert, and D. Baeriswyl, Z. 

Phys. B: Condens. Matter 103 (1997) 335. 
[19] S. Sorella and E. Tos atti, Europhys. Let t. 19 (1992) 699. 
[20] A. Vaezi, X. G. Wen. larXiv: 1010.57441 (2010). 
[21] T. Egami, et al. Science, 768, 1305 (1994); M. Fabrizio, A. O. 

Gogolin, A. A. Nersesyan, Phys. Rev. Lett.83 2014 (1999); M. 

Hafez, S. A. Jafari, Sh. Adibi, F. Shahbazi, Phys. Rev. B, 81 

(2010) 245131. 

[22] K. Bouadim, N. Paris, F. Hebert, G. G. Batrouni, R. T. Scalettar, 
Phys. Rev. B 76„ 085112 (2007). 



5 



[23] L. Craco, P. Lombardo, R. Hayn, G. I. Japaridze, E. Muller- (2009) 4479. 

Hartmann, Phys. Rev. B. 78, 075121 (2008). [25] D. Haberer, et. al. Physica Status Solidi (b), 248 (2011) 2639. 

[24] M. Hafez, S. A. Jafari, M. R. Abolhassani, Phys. Lett. A 373 



